/* coup__final.do 
Rabinowitz and Jargowsky 2017
*/


	* Cox Regressions 
			
		quietly stcox rps7 ib3.colony, vce(cluster ccode)
		estimates store m1
		
		quietly stcox rps7 ib3.colony ib1.epoch ib0.leftist prevcoup, vce(cluster ccode)
		estimates store m2
			
		quietly stcox rps7 ib3.colony ib1.epoch ib0.leftist prevcoup polity2, vce(cluster ccode)
		estimates store m3
			
		quietly stcox rps7 ib3.colony ib1.epoch ib0.leftist prevcoup polity2 ethfrag, vce(cluster ccode)
		estimates store m4
		
		quietly stcox rps7 ib3.colony ib1.epoch ib0.leftist prevcoup polity2 ethfrag relpol lpcgdp urbgrow, vce(cluster ccode)
		estimates store m5
		
		quietly stcox rps7 ib3.colony ib1.epoch ib0.leftist prevcoup polity2 ethfrag relpol lpcgdp urbgrow minpct oilpct, vce(cluster ccode)
		estimates store m6
		
		quietly stcox rps7 ib1.epoch ib0.leftist prevcoup polity2 i.ccode, vce(cluster ccode)
		estimates store m7
		
		quietly stcox rps7 ib1.epoch ib0.leftist prevcoup polity2 lpcgdp urbgrow minpct oilpct i.ccode, vce(cluster ccode)
		estimates store m8
		
	
		
		* Table 2
		esttab m1 m2 m3 m4, mtitle eform ///
			stats(N N_sub N_clust r2_p ll_0 ll chi2 , ///
			lab("Regime-Years" "Regimes" "Countries" "Pseudo R2" ///
			"Inital LL" "Final LL" "Chi-Squared")) ///
			star(* 0.05 ** 0.01 * 0.001) order(rps7) ///
			label drop(3.colony 1.epoch 0.leftist) ///	
			addnote("Base is French Colony, Democracy, Independence Period (pre-1980)")

		* Table 3
		esttab m5 m6 m7 m8, mtitle eform ///
			stats(N N_sub N_clust r2_p ll_0 ll chi2 , ///
			lab("Regime-Years" "Regimes" "Countries" "Pseudo R2" ///
			"Inital LL" "Final LL" "Chi-Squared")) ///
			star(* 0.05 ** 0.01 *** 0.001) order(rps7) ///
			label drop(3.colony 1.epoch 0.leftist *.ccode) ///	
			addnote("Base is French Colony, Democracy, Independence Period (pre-1980)")
		
		
		
//////////////////////////////////////////////////////////////////// END OF DO FILE///////////////////////////////////////////////////////


